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CLASFYT PROGRAM 
1.0 SCOPE 

This document contains the description of the CLASFYT program. The purpose 
of the program is as follows: 

To use a training field to determine trajectory coefficients for each 
band of radiance values. Five acquisitions of radiance values are 
used to create a composite set of tolerance values for each of the 
4 bands. 

To use these trajectory coefficients to classify each pixel in a test 
field to be the same agricultural product, not the same product or a 
screened pixel. 


2.0 APPLICABLE DOCUMENT'i 


The following documents form a part of this specification: 

AD 63-2457-3308-01 Transferring of Badhwar Software. 

JSC-16274, LEC-14064 Technical Memorandum Implementation of Badhwar 
Classification of Corn/Soybean Segments, G. D. Badhwar, W. W. Austin 
and J. C. Carnes, Agristars Report SR-jO-04001, JSC 16829. 

AD NAS 9-15200 Technical Memorandum Format Specification for LACIE 
(Phase III) and Accuracy Assessment Computer Data Products. 

M. J. Duggin, 9th Symposium on Remote Sensing of Environment, ERIM, 
Michigan (1974). 
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3.0 SYSTEM DESCRIPTION 


3*1 hardware description 

The software for CLASFYT Is operational on the IBM 3031 computer at Purdue, 


3-2 CLASFYT PROCESSOR SYSTEM FLOWCHART 

The system level data flow diagram Is shown In figure 3.1. 


The following system flowchart assumes the use of the EXEC file described 
■fn Appendix C. 


Aciiuisition data I'ilus (5) 
on multifile tape(s) 



(Segment number) 


Summary 
report on 
I disk and_^ 

LEciotef' 


.isrr A 


(Seg ment number) CLA SS A 


^ , Qassified 
-BH pixel file 




Figure 3.1 




CLASFYT Hierarchy Diagram 


CLASFYT (CLASFYT)*— CROPTM (CROPTM)i 


-MOVSYM 
1 — NUMBFL 
— NUMBR 


(CROPTM) 

(CROPTM) 

(CROPTM) 


p-RDEXCT (RDEXCT) 

— RDFLES (RDFLES) — RDVERT (RDFLES) 
— RECT (RECT) 

SDEV (SDEV) 

— -TRJFIT (TRJFIT)i— -F2ERIV (F2ERIV) 

— F2NCTN (F2ERIV) 
— FCHISQ (F2ERIV) 
— MATINV (F2ERIV) 
■-CHNFIT (CHNFIT)— FGERIV (F2ERIV) 

— FGNCTN (F2ERIV) 
■—FCHISQ (F2ERIV) 
--MATINV (F2ERIV) 


*Parenthese enclose the name of the file which contains the source statements 
for the named program which precedes them: PROGRAM NAME (SOURCE STATEMENT 
FILE NAME). 
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Following Is a brief description of CLASFYT program Inputs and outputs. 
For a more detailed description refer to section 3.4. 


Acquisition data files 


User Information file 


Data from RT&E data base 


Summary report 


Classified pixel file 


- Header record and 4 channels of data for 
eac^’ 5 acquisitions (See Section 3.4.3). 

- User Information file supplying crop, channel 
processing order. Initial guesses for each 
of 4 channels, training field vertices and 
optional test field vertices (See Section 
3.4.1). 


file - Segment number, acquisition dates and 
Landsat numbers (See Section 3.4.2). 

- Report of date and time of processing, crop, 
vertices for training field, acquisition 
dates, all statistical data created by 
CLASFYT and time taken to execute the 
program. 

- A file consisting of a header record and 
117 data records. Each data record 
contains the classification code for each 
of the 196 pixels In a line of acquisition 
data (See Section 3.4.4). 


3.3 SORVIARE DESCRIPTION 

This program classifies each pixel in a segment into one of several classes. 
The statistical characteristics of the classes are defined using the pixels 
in training fields; there is one training field for each clas« up to a maximum 
of ten classes. If a pixel is not classified as any of the input classes it 
is put into a special class called "unclassified". For a description of the 
method used see reference 3 in Section 2.0. 

For each training field the program does the following; 

(1) Chooses the first 40 pixels in the training field and for 
these pixels calculates the mean and standard deviation for 
each channel of each acquisition. (The other pixels in tho 
training field are not used.) 

(2) Eliminates any of the 40 pixels that have values (landsat 
counts) more than thr.i-e times the standard deviation away 
from the mean in any channel of any acquisition. 

(3) Recalculates the mean and standard deviation for the remaining 
pixels . 

(4) Calls the routine TRJFIT which estimates the values of A, a, 3, 
and t„ which give the best fit of the function 

Ftt) =a(|)“ 0 '^“'■‘-*0 

and ° 

F(t) = A for t < t^ 

to the time trajectory of the means in a particular channel. 

This is done for each channel. The fitting routine is based 
on the statistics calculated in step 3 above and a set of 
user-defined initial guesses for the values of the parameters 
A, a, 6 and t^ (Section 3.1). 



3-4 


c *»r« 


1[SJ Estimates the standard deviation due to sensor and scene noise 
(Ref. 3, Section 2.0) for each channel In each acquisition. 

(6) Calls the routine CHNFIT for each channel. This routine fits 

F(t) to the time trajectory of each of the remaining pixels 

of the original set of 40. In doing tills the values of A, 

a and d obtained In step 4 above are used, but a new value 

2 

of tg for each pixel Is determined. Also a value for x 
which describes the goodness of fit Is calculated for each 
pixel. This calculation Is based on estimates of the 
"measurement noise" for each acquisition. This quantity Is 
taken to be the larger of the standard deviations calculated 
In steps 3 and 5. 

(7) The means (over pixels) and standard deviations of the’ x^, values 
are calculated for each channel and used to define a x" 
threshold value for each channel. 

. (8) The means and standard deviations of the t^ values 
"afe calculated for eac^ and used to define “ 

an upper and lower t^ threshold value for each channel. 

Next, each pixel In the segment Is classified using the results obtained In 
steps 1 to 8 above. For each pixel the following procedure Is followed: 

0) The pixel Is tested to determine If It should be classified as 

class 1. To do this each channel of data Is processed separately 

In the order specified by the user. The default channel order 

Is 2, 3, 4, 1. For each channel the routine CHNFIT Is called, 

2 

As described above, this routine provides an estimate of x 

2 

and t^ for this channel. If the value of x Is above the 
threshold values for t^ calculated In step 8 above, then the 
pixel Is not cl a«1 f ied“a‘sTTass' 1 . If"the”p1xeT"Ts 
not rejected In any of the 4 channels then It Is classified 
as class 1 . 
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(2) If the pixel Is rejected as class 1 In any channel, then It Is 
tested in the same way to determine if it belongs to class 2. 

If it is rejected as class 2 then it is tested for class 3, etc. 
If it is rejected for all classes the it is put in a special 
class called "unclassified*^. 
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3.4 FILE DESCRIPTION 


3.4.1 USER INFORMATION FILE (.Segment number) USER! A 

The following records describe the order of channel processing, the name, 
initial guesses for the 4 channels and the symbol for each of the crops, 
the training and test fields and allows the user to Input any number of 
free-form conments. 

In record types 1-9 the ''keyword'' Is In columns 1-4 and any parameter data 
is in columns 11 through 72. Numbers In a series are separated by commas; 
blanks are optional. 

Record types 10 and 11 which describe the fields are blank In columns 1 
through 10 and the field definition Is In columns 11-72. 

1. "COMMENT" record 

EXAMPLE: COMMENT THIS IS A TEST EXECUTION OF THE PROGRAM 

These records allow the user to Input any number of free-form comments. 

2. "CHANNEL" record 

EXAMPLE: CHANNEL 2, 3, 4, 1 

This card specifies the order In which the channels are to be evaluated. 
The default is 2, 3, 4, 1. 

3. " CROP " record (s) 

EXAMPLE: CROP CORN, S OYBEANS 

This record specifies the 4 character names of the agricultural products. 
The names begin with the first non-blank character and end at the fourth 
character or the comma separating the names. Any additional characters 
are ignored. 


4- “ Q IN ITIAL GUESS" record 

EXAMPLE: Q INIT 3,45, -3.65, -0.32, 1.50 

Q INIT 3.5, -1.0, -1.0, 1.2 

These records, specify the A, Alpha, Beta and t^ Initial guesses for 
channel 1 for each of the crops to be evaluated. 
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'5. "X INITIAL GUESS" record(s) 

EXAMPLE: X INIT 3.41, -9.61, -0.91, 1.50 (Corn) 

X INIT 3.5, -5.0, -1.0, 1.20 (Soybeans) 

These records specify the A, Alpha, Beta and t^ initial guesses for 
channel 2 for each of the crops to be evaluated. 

”Y INITIAL GUESS" record(s) 

EXAMPLE: Y INIT 3.08, 8.65, 1.01, 1.47 

Y INIT 3.5, 7.0, 1.0, 1.2 

These records specify the A, Alpha, Beta and t^ Initial guesses for 
channel 3 for each of the crops to be evaluated. 

7. "Z INITIAL GUESS" record(s) 

EXAMPLE: Z INIT 2.91, 11.17, 1.3, 1.05 

Z INIT 3.5, 7.0, 1.0, 1.2 

These records specify the A, Alpha, Beta and tg initial guesses for 
channel 4 for each of the crops to be evaluated. 

8. "SYMBOL" record 
EXAMPLE: SYMBOL 195, 226 

This record specifies the symbols for the crops in the same order they 
are listed in the CROP record. Each symbol is defined by its base 10 
number. These numbers are designated by crop codes produced by the data 
management section at J.S.C. and vary according to crop and year. 

9 . "*END" record 
EXAMPLE: *END 

This record specifies the end of the user supplied data except for the 
definition of the training field and the test field. 

Record type 10 must follow the *END record. 

Record type 11 is optional, the default is the whole scene. If record 
type 11 is present it must follow record type 10. 
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10. "TRAINING FIELD" record 


(The training field record must follow the "*END" record and proceed 
the test field record. 

EXAMPLE; (1,1), (97,90 ), (125, 90), (100,100), (80,100) 

This record specifies the coordinates of the corners of the training 
field. The format of the line Is as follows: 

(1,1) Duminy variable to be used for sample skip factor and 
line skip factor; not currently used. 

(Xi,Yi) Upper leftmost corner of training field expressed as 
(sample, line). 

Upper rightmost corner of training field expressed as 
(sample, line). 

(Xa.Y^) Lower rightmost corner of training field expressed as 
(sample, line). 

(^ 4 »Y/^) Lower leftmost corner of training field expressed as 
(sample, line). 

11. "TEST FIELD" record 

(The test field record. If specified, must follow the training field 
record. The default Is the whole scene. 

EXAMPLE: (1,1), (1,1), (196,1), (196,117), (1,117) 

This optional record specifies the coordinates of the corners of the 
(test field. The default is the whole segment. The format Is the same 
as the format for the training field record. 
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3.4.2 ACQUISITION DATA FILES 

The following records of thTs' acquisi ti on s data' fi 1 es ’are* friput to" thi" ~~ 
CLASFYT program and are read by RDFLES. 

These files each contain 4 channels of LANDSAT data for one acquisition. 

If an acquisition is used twice the corresponding file will be duplicated. 
Note that these files are in universal format. Record 1 is read with 
a 17(180A1) format statement; records 2-118 are read with a 5(180A1) format 
statement. 


Record 

Bytes 

Contents 

1 

1 -3060 

Header record in universal format. 

2 

1-72 

Filler. 


73-268 

Gray level values for the 196 samples of band 1 
for line 1 . 


296-464 

Gray level values for the 196 samples of band 2 
for 1 ine 1 . 


465-660 

Gray level values for the 196 samples of band 3 
for line 1 . 


661-856 

Gray level values for the 196 samples of band 4 
for 1 ine 1 . 

3 

1-856 

Same data format as record 2; gray level values 

• 


for bands 1-4 for line 2. 

• 

• 

N 

1-856 

Same data format as record 2; gray level values 

• 


for bands 1-4 for line N-1. 

• 

• 

118 

1-856 

Same data fromat as record 2; gray level val.ues 


for bands 1-4 for line if?. 
^ « 
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3.4.3 DATA FROM RT&E DATA BASE FILE 


LANOSN FILE D 


The following records describe the segment number and the date and Landsat 


number 

for each of the 5 acquisitions. 


Record 

Bytes 

Format 

Contents 

1 

1-4 

A4 

Segment number. 

2 

1-35 

5(2X^5) 

Acquisition dates. 

3 

1-20 

514 

Landsat numbers. 
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3.4.4 CLASSIFICATION FILE (Segment number) CLASS A 

The following records are output from the CLASFYT program. 
Format 


Record 

Bytes 

Format 

Contents 

1 

1-3060 

17(180A1) 

Header record in Universal format. 

2-118 

1-72 

17(180A1) 

72 characters of blank filler. 


73-268 


196 characters representing 1 line 
of classified data. 

Record 

Bytes 

Contents 


1 

1-3060 

Header records 

in Universal format. 

2 

1-72 

Blank filler. 



73 

Classification 

symbol for pixel 1 in line 1. 



N+72 

• 

Classification symbol for pixel N in line 1. 


• 

268 

Classification symbol for pixel 196 in line 1. 

3 

1-268 

Same data format as record 2; 72 characters of blank 
filler and classification symbols for 196 pixels in 

• 


line 2. 

• 

• 

N+1 

1-268 

Same data format as record 2; 72 characters of blank 
filler and classification symbols for 196 pixels in 
line N. 
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Classification file (continued). 


Record 

118 


Byte Contents 

1-268 Same data format as record 2; 72 characters of blank 

filler and classification symbols for 196 pixels In 
line 117, 
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3.5 DETAILED SOFTWARE DESCRIPTION 

3.5.1 CLASFYT 
Purpose 

CLASFYT uses a user determined training field In a full scene to determine 
trajectory coefficients for each chanM. of radiance value s T " These coefficients" 
are used to classify each pixel In a test field to be a specific one of the 
crops of Interest or unclassified. 

Linkages 

CLASFYT calls CROPTM, RDEXCT, CPTIME, RECT, SDEV, TRJFIT, CHNFIT and RDFLES. 

Interface 
Calling sequence: 

Not applicable (an EXEC that can be used to call CLASFYT Is described 
In Section 4.0) . 

Calling sequence parameters: 

Not applicable. 

■Labeled COMMON parameters: 


Full description of labeled 

COMMON blocks 
Element 

are contained in Appendix A. 
Input/ 

Label 

Variable 

Position 

Output 

Description 

/LINES/ 

LINENO 

2 

0 

Line number. 

/LINES/ 

IDATA(856,8) 

3 

I 

One line of ^ 
acquisition data per 
file. 

/LINES/ 

NOAcqs 

4 

0 

Cumber of acquisitions. 

/USER/ 

ICROP(IO) 

1 

I 

Names of crops. 


Label 

/USER/ 

/USER/ 

/USER/ 

/USER/ 

/USER/ 

mm~ 

/USER/ 

/USER/ 

/MISC/ 

/MISC/ 

/MISC/ 

/MISC/ 


Blank COMMON 
None. 
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Variable 

Position 

Input/ 

Output 

Description 

ZQ(4) 

2 ' 

I 

Arrays of Coefficients 

ZX(4) 

3 

I 

A Alpha, Beta and t 
for channel 1 . ° 

Array of coefficients 

ZY(4) 

4 

I 

A, Alpha, Beta and t 
for channel 2. ° 

Arrays of coefficients 

ZZ(4) 

5 

I 

A, Alpha, Beta and t 
for channel 3. ° 

Arrays of coefficients 

NMAX 

6 

I 

A, Alpha, Beta and t 
for channel 4. ° 

Maximum number of 

70RDRCC4)^ -- 
NOCROP 

8 

I_„ 

I 

points In field 
(default Is 40). 

-Array-specl-fylng order 
channels are to be 
processed. 

Number of crops. 

ISYMBL(IO) 

9 

I 

Symbols for crops. 

NOACq 

1 

I 

Number of acquisitions. 

ACQDTC5) 

2 

I 

Acquisition dates. 

KSE6M 

3 

I 

Segment number. 

LANDST(5) 

t 

4 

I 

Landsat numbers. 


parameters: 


Inputs 


Unit 

Type 

Description 


7 

Seq. data 

Data from RT&E data base (See 

Section 3.4.3), 

11 

Seq. data 

Acquisition data (See Section 

3.4.2). 

12 

Seq. data 

Acquisition data (See Section 

3.4.2). 

13 

Seq. data 

Acquisition data (See Section 

3.4.2). 

14 

Seq. data 

Acquisition data (See Section 

3.4.2). 

15 

Seq. data 

Acquisition data (See Section 

3.4.2). 

21 

Seq. data 

User Information (parameters) 

(See Section 3.4.1). 

30 

Seq. data 

Scratch file to reformat data. 

Outputs 

Unit 

Type 

Description 


6 

Report 

Summary report. 


8 

Seq. data 

Classified pixel file. 


30 

Seq. data 

Scratch file to reformat data (See Section). 


Storage requirement 
Not applicable. 

Description 

CLASFYT calls CROPTM to read the user specified parameters crops of interest, 
initial guesses for coefficients A, Alpha, Beta, t^, and symbols for each 
crop. 

CLASFYT calls RDFLES to read the header record and the parameters for the 
training field. 
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CLASFYT calls RDHLES to read the header record and the parameters for the 
training field. 

CLASFYT then calls CPTIME to get the current time and date and writes the 
report heading. 

CLASFYT calculates the minimum and maximum boundaries for pixels In the test 
field, calculates days elapsed from the beginning of the year as a multiple 
of 100 for each acquisition date and sets a switch If an acquisition date Is 
used twice. 

The rectangular field bounded by the minimum and maximum pixel numbers and 
line nu oers is then processed using the following steps: 

A line of radiance values Is read. 

Each pixel In the field is checked to determine If It Is In the training 
field. 

The^alues for the pixels In the trai ning fields are saved In a matrix . 

When sufficient training field pixels (40) have been Identified/accumulated, 
the mean and standard deviation for each of the four channels of radiance 
values are calculated for the selected training field pixels. 

The radiance value for each of the four channels for the selected training 
field pixels Is compared with the corresponding computed channel mean and 
channel standard deviation. Any selected training field pixel whose radiance 
value for any channel differs from the corresponding channel mean be more 
than 3 channel standard deviations Is eliminated from the training field 
matrix. 

The mean and standard deviation are re-calculated for the new group of 
training field pixels and are printed on the report. 
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CLASFYT calls TRJFIT to calculate CHI-square coefficients for each channel. 
The final coefficients are calculated using the CHI-square values and the 
constants from Duggins' paper (1974) for scene noise and CHNFIT is called 
for each channel. The initial guesses and final coefficient are written 
on the report for each channel. 

The header record is constructed using input from the program constants » 
the user parameter file and the current date. The header record is wifitten 

as the first record on the new classification file, 

^ .. 

CLASFYT then evaluates the test field using the following steps: 

A line of pixel data is read. 

CHNFIT is called for each of the pixels for channel 4, channel 2, 
channel 3 and then channel 1 unless a different sequence is specified 
by the user. 

The pixel is classified as a different product: if the CHI square 
value calculated by CHHFIT does not fall within the limits calculated 
for the corresponding channel. 

The line of classified pixels is written to the classification file. 

The threshold values and the counts of values classified are written to the 
report. 

FI owchart 
Not applicable. 


Listing 

See Appendix D for program. 


3,5.2 SUBROUTINE CROPTM 
Purpose 


CROPTM reads and analyzes cards describing the order of channel processing 
and the name of the crops of Interest, Initial values for A, Alpha, Beta and 
t^, and symbols for each crop. 

LI nkaqes 

CROPTM Is called by CLASFYT, CROPTM calls MOVSYM, NUMBFL, and NUMBR. 

Interface 

Calling sequence: 

CALL CROPTM (lERROR). 

Calling Sequence parameters: 


Input/ 


Arqument Output 

Description 

I ERROR 0 

0 - No error; 

Function value: 

Nolf applicable. 

1 - error. 


Labeled COMMON parameters; 


Full description of labeled COMMON blocks are 

Element Input/ 

contained In Appendix A 

Label 

Variable 

Position 

Output 

Description 

/USER/ 

ICROP(IO) 

1 

0 

Names of crops. 

/USER/ 

ZQ(4,10) 

2 

0 

Arrays of coefficients 
A, Alpha, Beta and t 
for channel 1 . ° 

/USER/ 

ZX(4,10) 

3 

0 

Arrays of coefficients 
A, Alpha, Beta and for 
channel 2. 

/USER/ 

/ 

ZY(4,10) 

4 

0 

Arrays of coefficients 
A, Alpha, Beta, and t 
for channel 3. ° 
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/USEFV 

QZ(4,10) 

5 

j3 Arrays of coefficients 

A, Alpha, Beta and t 
for channel 4. 

/USER/ 

MAXPTS 

6 

0 Maximum points. 

/USER/ 

I0R0RCC4) 

7 

0 Array describino order 

channels are to be processed. 

/USEIV 

NOCROP 

8 

0 Number of crops . 

/USER/ 

ISYMBLCIO) 

9 

0 Array of symbols. 


Blank COf^ON parameters: 

None. 

Inputs 

Unit Type Description 

21 Seq. data User information file (See Section 3.4.1), 

Outputs 

None. 

Storage requirement 
Not applicable. 

Description 

Each record in the user defined parameter file is read and checked to have a 
valid control word in characters 1 through 4 defining the type of information , 
in characters 11 through 80. If the control word in invalid, an error 
message is written; otherwise, the information is saved in the appropriate 
comnon block element. 


FI owcha rt 
Not applicable. 

Listing 


See Appendix D for program. 




3.5.3 FUNCTION MOVSYM 
Purpose 

MOVSYM moves characters from array CARD to a series of left-justified 
words in array ITEMP. 

Li nkaqes 

MOVSYM is called by SETUPT. 

Interface 

Calling sequence: 

KCUNT = MOVSYM (CARD, ITEMP). 

Calling sequence parameters: 

Input/ 

Argument Output 
CARD I 

ITEMP 0 

Function value: 

Name Description 

KOUNT Count of symbols. 

Labeled COMMON parameters: 

None. 

Blank COMMON parameters: 

None. 

Inputs 


Description 

80 characters of user supplied data in 
character format. 

Array of left justified words or symbols. 


None. 


Outputs 

None. 


Storage requirement 
Not applicable. 

Description 

MOVSYM separates the characters in array CARD starting at character 11 and 
creates a series of left-justified words in array ITEMP. 

Flowchart 
Not applicable. 


Listing 

See Appendix D for program. 


3.5.4 FUNCTION NUMBR 
Purpose 


ORIGINAL PAGE IS 
OF POOR QUALITY 


NUMBR changes integers from character array CARD to computational numbers in 
array NUMVEC. 

Linkages 

NUMBR is called by CROPTM. 

Interface 

Calling sequence: . 

KOUNT = NUMBR (CARD, NUMVECyT~~ 


Calling sequence parameters: 


Argument 

Input/ 

Output 

Description 

CARD 

I 

80 characters of user supplied data in character 
format. 

NUMVEC 

0 

Array of integer computational numbers. 

Function value 
Name 

• 

• 

Description 



KOUNT Number of computational integers. 

Labeled COMMON parameters: 

None. 

Blank COMMON parameters: 

None. 

Inputs 

None. 

Outputs 

None. 
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storage requirement 
Not applicable. 

Description 

NUMB separates the characters In array card starting at the 11th character 
and creates a series of computational Integers. 

FI owchart 
Not applicable. 

Listing 

See Appendix D for program. 
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3.5.5 FUNCTION NUMBFL 
Purpose 

NUMBFL changes floating point numbers from characters in array CARD to 
computational numbers in array XNMVEC. 

Li nkaqes 

NUMBFL is called by CROPTM,. 

Interface 

Calling sequence: 

KOUNT = NUMBFL (CARD, XNMVEC). 

Calling sequence parameters: 

Input/ 

Argument Output 
CARD I 

XNMVEC .0 

Function value: 

Name Description 

KOUNT Count of floating point numbers. 

Labeled COMMON parameters: 

None. 

Blank COMMON parameters: 

None . 

Inputs 


Description 

80 characters of user supplied data in • 
character fonnat. 

Array of floating point computational 
number. 


None. 


Outputs 

None. 


Storage requirement 
Not applicable. 

Description 

NUMBFL evaluates the characters in array CARD starting at character 11 and 
creates a series of computational floating point numbers. 

FI owchart 
Not applicable. 


Listing 

See Appendix D for program. 


3.5.6 SUBROUTINE RDEXCT 
Purpose 


RDEXCT reads a file composed of the acquisition information from the LARS 
data base. 

Linkages 

RDEXCT is called bv CLASFYT. 

Interface 

Calling sequence: 

CALL RDEXCT. 

Calling sequence parameters: 

None. 


Function value: 

Not applicable. 


Labeled COMMON parameters: 

Full description of labeled COMMON blocks are contained in Appendix A 


Label 
/MI SC/ 
/MISC/ 
/MISC/ 
/MISC/ 


Variable 

NOACQ 

ACQDT(5) 

KSEGM 

LANDST(5) 


Element 

Position 

1 

2 

3 

4 


Input/ 

Output 

I 

I 

I 

I 


Description 
Number of acquisition 
Acquisition dates. 
Segment number. 
Landsat numbers. 


Blank COMMON parameters; 
None. 
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Inputs 


Unit Type Description 

7 Seq. data Acquisition information from RT&E data base. 

(See Section 3.4.3). 

Outputs 

None. 

Storage requirement 
Not applicable. 

Description 

RDEXCT reads a fixed format file defining the segment number, acquisition 
dates and Landsat numbers. RDEXCT sets the number of acquisitions to 5. 

Flowchart 
Not applicable, 

Li s ti ng 

See Appendix D for program. 
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3.5.7 SUBROUT INE RDFLES 
Purpose 

RDFLES reads 1) the header record and one set of vertices or 2) one line of 
radiance values. 

Linkages 

READFL is called by CLASFYT. 

Interface ~ ~ — - — 

Calling sequence: 

CALL RDFLES (ITYPE, HDR, XFD, XFD, lERROR) . 

Calling sequence parameters: 

Input/ 

Argument Output 

ITYPE I 

HDR 0 

XFD(4) 0 

YFD(4) 0 

lERROR 0 

Function value: 

Noit applicable. 

— — - ~ - j 


Descri ption 

1 to read header* from each acquisition file 
and field definition; 2 to read 1 line of 
radiance values from each acquisition file. 

Header record buffer. 

Array of sample vertices defining field. 

Array of sample lines defining field. 

0 if no error; 

1 if error. 
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Labeled COMMON parameters: 

Full description of labeled COMMON blocks are contained In Appendix A. 


Label 

Variable 

Element 

Position 

Input/ 

Output 

Description 

/LINES/ 

NXTLNE 

1 

0 

Pointer to next line. 

/LINES/ 

LINENO 

2 

I 

Number of line to be 
read. 

/LINES/ 

IDATA(856,8) ,3 

0 

Line of radiance 
values for each 
acquisitions. 

/LINES/ 

NOACQS 

4 

0 

Number of acquisitions 

Blank COMMON parameters: 




None. 





Inputs 





Unit 

Type 

Description 



11 

Seq. data 

Acquisition data 

for date 

1 (See Section 4.3.2). 

12 

Seq. data 

Acquisition data 

for date 

2 (See Section 4.3.2). 

13 

Seq. data 

Acquisition data 

for date 

3 (See Section 4.3.2). 

14 

Seq. data 

Acquisition data 

for date 

4 (See Section 4.3.2). 

15 

Seq. data 

Acquisition data 

for date 

5 (See Section 4.3.2). 

21 

Seq. data 

Field definition. 


Ouputs 

None. 




f 


Storage ‘requirement 
Not applicable. 


3-31 


Description 

REAOFL reads the header record from the universal formatted tape and the 
vertices for one field if ITYPE«*1 . If ITYPE°2, READFL reads one line of 
radiance values. 

FI owcharts 
Not applicable. 

Listing 

See Appendix D for program. 
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3.5.8 SUEkOUTINE RDVERT 
Purpose 

RDVERT reads one set of vertices. 


Linkages 

RDVERT Is called by CLASFYT. 

Interface 

Calling sequence; 

CALL RDVERT (XFLD, YFLD, lERROR). 


Calling sequence parameters: 

Input/ 

Argument Output 
XFLD(4) 0 

YFLD(4) 0 

lERROR 0 


Description 

Array of sample vertices defining field. 
Array of line vertices defining field. 

0 - no errors; 

1 - error. 


Function value: 

Not applicable. 


Labeled COMMON parameters: 
None. 


Blank COMMON parameters: 
None. 
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Inputs 

Unit Type Description 

21 Seq. data User specified data file (See Section 4.3.1). 

Outputs 

None. 

Storaqe requirement 
Not applicable. 

Description 

RDVERT reads the vertices for one field. 

Flowchart 
Not applicable. 

Listing 

See Appendix D. 
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3.5.9 SUBROUTINE RECT 
Purpose 

RECT determines If a pixel lies within the specified vertices of a specific 
area. 


Li nkaqes 

RECT is called by CLASFYT. 

Interface 
Calling sequence: 

CALL RbCT CLINE, PIXEL, YY, XX, ICOR) . 


Calling sequence parameters: 

Input/ 

Argument Output Description 


LINE I 

PIXEL I 

YY(4) I 

XX(4) I 

ICOR 0 


Current line number of pixel. 

Current sample number of pixel. 

Array of sample vertices defining field. 
Array of line vertices defining field. 

0 if pixel in field; 

1 if pixel not in field. 


Function value: 

'Not applicable. 


Labeled COMMON parameters: 
None. 

Blank COMMON parameters: 
None. 










Inputs 

None. 

Outputs 

None. 

Storage requirement 
Not applicable. 

Description 

The pixel is checked to be a boundary pixel and considered in the training 
field if it is a boundary pixel. If the pixel is within the minimum and 
maximum limits of tl i specified area, the angles created between lines 
drawn between the pixel and boundary pixel are calculated and summed. If 
the sum of the angles is 360; the pixel is considered to be in the specified 
field. 


Flowchart 
Not applicable. 

Listing 


See Appendix D for program. 


3.5.10 SUBROUTINE SDEV 
Purpose 

SDEV calculates the mean and standard deviation of a vector. 


Linkages 

SDEV is called by CLASFYT. 

Interface 

Calling sequence: 

CALL SDEV CX, N, XB, S). 


Calling sequence parameters: 

Input/ 

Argument Output Description 


X(N) 

N 

XB 

S 


I Array of values. 

I Number of values in array X. 

0 Mean of values in array X. 

0 Standard deviation of values in array X 


Function value: 

Not applicable. 


Labeled COMMON parameters: 
None. 


Blank COMMON parameters: 
None. 


Inputs 

None. 
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storage Requirement 
Not applicable. 


Description 

SDEV calculates the mean (XB) and standard deviation (S) of vector X with 
N points. 

Flowchart 
Not applicable. 

Listing 


See Appendix D for program. 
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3.5.11 SUBROUTINE TRJFIT 
Purpose 

TRJFIT makes a least-square fit to a non-linear function to calculate 
coefficients for the curve. 


Linkages 

TRJFIT is called by CLASFYT. 

Interface 
CalMng sequence: 

CALL TRJFIT (X, Y, SIGMAY, NFILES, NTERMS, MODE, ESTIM, DELTAA, SIGMAA, 
FLAMDA, YFIT, CHISQR) . 


Calling sequence parameters: 

Input/ 

Argument Output 


X I 

Y I 

SIGMAY I 

NFILES I 

NTERMS I 

MODE I 

ESTIM 1,0 

DELTAA I 

SIGMAA 0 

FLAMDA 0 

YFIT 0 

CHISQR 0 


Description 

Array of Julian day/100 for acquisitions. 
Array of dependent variables. 

Array of errors on Y. 

Number of acquisitions. 

Number of degrees of freedom. 

Not used, was switch for calculation weight. 
Array of coefficients A, Alpha, Beta and t^. 
Array of increments for A. 

Standard deviations . 

Proportion of gradient search. 

Array of fitted values. 

Residual CHI-square value. 




Function value: 

Not applicable. 

Labeled COMMON parameters 
None. 

Blank COMMON parameters: 
None. 

Inputs 

None. 

Outputs 

None. 

Storage requirement 


Not applicable. 


Description 


TRJFIT calculates alpha and beta matrices, using the least-square method, 
for defining all parameters of the curve for the selected pixels in the 
training field. 

The CHI square value of the starting point is evaluated and the modified 
curvature matrix is inverted to find the new parameters. 

The convergence depends on the CHI-square value and on the size of step 
FLAMDA which controls the location of the emergence day. 

If the CHI-square value increases, FLAMDA is increased and the process 
is re-initiated. If the CHI-square value decreases, the parameters and 
uncertainties are evaluated. 

Flowchart 
Not applicable. 

Listing 

See Appendix D for program. 
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3.5.12 SUBROUTINE F2ERIV 
Purpose 

F2ERIV computes the derivative of the function of each of the X values. 
Linkages 

F2ERIV Is called by TRJFIT. 

Interface 
Calling sequence: 

CALL F2ERIV (X, I, ESTIM, DELTAA, NTERMS, DERIV) . 


Calling sequence parameters: 

Input/ 

Argument Oubut 


X I 

I I 

ESTIM 1,0 

DELTAA I 

NTERMS I 

DERIV 0 


Description 

Array of time points for acquisition. 

ft 

Current Index In calling program. 

Array of coefficients A, Alpha, Beta and t^. 

An ’ay of Increments for A. 

Number of degrees of freedom. 

Array of derivatives with respect to parameters 


ESTIM. 


Function value: 

Not applicable. 


Labeled COMMON parameters: 
None. 
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Blank COMMON parameters: 
None. 

Inputs 

None. 

Outputs 

None. 

Storage requirement 
Not applicable. 


Description 

The radiance value and t^ are set positive. Y is calculated as the radiance 
value/t^. If Y is less than or equal to 1, the 4 derivatives are set to 
very small numbers; otherwise, the derivatives are set to calculated values. 
If the CHI-square value decreases, the process is reinitiated until 
convergence is obtained; the parameters and uncertainties are evaluated. 

Flowchart 
Not applicable. 

Listing 


See Appendix D for program. 


3.5.13 FUNCTION F2NCTN 
Purpose 

F2NCTN evaluates a function of alpha, beta, to and time. 
Linkages 

F2NCTN Is called by TRJFIT. 

Interface 

Calling sequence: 

APPROX = F2NCTN (X, I, COEFS). 


Calling sequence parameters: 

Input/ 

Argument Output 
X I 

I I 

COEFS I 


Description 

Time array In days elapsed In year/ 100. 
Index to value In X array. 

Array of estimates A, Alpha, Beta and t^. 


Function value: 

Name Description 

APPROX Evaluation of function of Alpha, Beta, t^ and X. 


Labeled COMMON parameters: 
None. 


Blank COMMON parameters: 
None. 




Inputs 

None. 


Outputs 

None. 
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storage Requirement 
Not applicable. 


Description 

F2NCTN sets the variable to positive and then evaluates the following 
equation Y « radiance value/titne. If Y Is less than or equal to 0, 

F2NCTN = radiance value, else Y = radiance value + alpha + alog (Y) + Beta * 
To^ * (l-t1me^). 

Flowchart 
Not applicable. 

Listing 

See Appendix D for program. 
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3.5.14 FUNCTION FCHISQ 
Purpose 

FCHISQ calculates the CHI square value per degree of freedom. 
Linkages 

FCHISQ Is called by TRJFIT and CHNFIT. 

Inter'*^ace 
Calling sequence: 

GOOOFT = FCHISQ (Y, SIGMAY, NPTS, NFREE, MODE, YFIT, WEIGHT). 
Calling sequence parameters: 


Input/ 


Argument 

Output 

Description 

Y(NFILES) 

I 

Array of dependent variables. 

SIGMAY(NFILES) 

I 

Array of errors on Y. 

NPTS 

I 

Number of acquisitions. 

NFREE 

I 

Number of degrees of freedom. 

MODE 

I 

Not used. 

YFIT(NFILES) 

I 

Array of fitted values. 

WEIGHT(NFILES) 

I 

Array of weighting factors. 

Function value: 

Name Description 


GOODFT CHI-square 

1 value 

per degree of freedom. 


Labeled COMMON parameters; 
None. 
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fc 

f 

\ 

Blank COMMON parameters: 

None. 

Inputs 
None . 

I Outputs 

' None. 

I Storage Reg u remen t 

p Not applicable. 

Description 

FCHISQ calculates CHI square value = I (YFIT-Y)/SIGMA)^ 

Flowchart 
Not applicable. 

) Listing 

See Appendix D for program. 
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3.5.15 SUBROUTINE MATINV 
Purpose 

MATINV inverts a matrix of order NORDER. 
Li nkaqes 

MATINV is called by TRJFIT and CHNFIT. 

Interface 

Calling sequence; 

CALL MATINV (ARRAY, NORDER, DET). 

Calling sequence parameters: 

Input/ 

Argument Output 

ARRAY (NORDER) 1,0 

NORDER I 

DET 0 

Function value: 

None. 

Labeled COMMON parameters: 

None. 

Blank COMMON parameters: 

None. 

Inputs 
None. 

Outputs 
None. 


Description 

Matrix. 

Order of matrix. 
Determinate of matrix. 
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storage requirement 
Not applicable. 


Description 

MATINV uses the Gauss-Jordon reduction techniques. 

Flowchart 
Not applicable. 

Listing 


See Appendix D for program. 


3.5 16 SUBROUTINE CHNFIT 
Purpose 

CHNFIT makes a least-square fit to a non-linear function to calculate 
coefficients for the curve. 

Linkages 

CHNFIT is called by CLASFYT. 

Interface. 

Calling sequence: 

CALL CHNFIT (X, Y, SI6MAY, NFILES, NTERMS, MODE, COEFS, DELTAA, SIGMAA, 
FLAMDA, YFIT, CHISQR, CON, CALP, CBET) . 

Calling sequence parameters: 


Argument 

Input/ 

Output 

Description 

X 

I 

Array of Julian day/lOO far acquisitions. 

Y 

I 

Array of dependent variables. 

SIGMAY 

I 

Array of errors on Y. 

NFILES 

I 

Number of acquisitions. 

NTERMS 

1 

Number if degrees of freedom. 

MODE 

I 

Not use,, was switch for calculation weight. 

COEFS 

1,0 

Array of coefficients A, Alpha, Beta and t^. 

DELTAA 

I 

Array of increments for A. 

SIGMAA 

0 

Standard deviations. 

FLAMDA 

0 

Proportion of gradient search. 

YFIT 

0 

Array of fitted values. 

CHISQR 

0 

Residual CHI-square value. 
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Input/ 

returned from TRJFIT (A) . 
returned from TRJFIT (Alpha), 
returned from TRJFIT (Beta) . 

Function value: 

Not applicable. 


Argument 

Output 

Description 

CON 

I 

Data point 1 

CALP 

I 

Data point 2 

CBET 

I 

Data point 3 


Labeled COMMON parameters: 

Full descriptions of labeled COMMON 


Element 

Label Variable Position 
/CHAN/ ACH 1 

/CHAN/ ALP 2 

/CHAN/ ABET 3 


Blank COMMON parameters: 
None. 

Inputs 

None. 

Outputs 

None. 


blocks are contained in Appendix A. 
Input/ 

Output Description 

0 A returned from TRJFIT. 

0 Alpha returned from TRJFIT. 

0 Beta returned from TRJFIT. 


Storage requirement 
Not applicable. 
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Description 

CHNFIT calculates alpha and beta matrices, using the least-square method, 
for defining all parameters of the curve for the selected pixels in the 
training field. 

The CHI square value of the starting point is evaluated and the modified 
curvature matrix is inverted to find the new parameters. 

The convergence depends on the CHI-square value and on the size of step 
FLAMDA which controls the location of the emergence day. 

If the CHI-square value increases, FLAMDA is increased and <-he process 
is re-initiated. If the CHI-square value decreases, the parameters and 
uncertainties are evaluated. 

Flowchart 
Not applicable. 

Listing 


See Appendix D for program. 


3.5.17 SUBROUTINE FGERIV 
Purpose 

FGERIV calculates the value of DERIV from the radiance value, alpha and 
beta. 


. Linkages 

FGERIV is called by CHNFIT. 


Interface 
Calling sequence: 

CALL FGERIV (X, I, COEFS, DELTAA, NTERMS, DERIV). 


Calling sequence parameters: 

Input/ 

Argument Output 


X I 

I I 

COEFS I 

DELTAA NA 

NTERMS NA 

DERIV 0 


Description 

Time array in days elapsed in year/100. 
Index to value in time array. 

Array of coefficients A, Alpha, Beta and t^ 
Not used. 

Not used. 

Value of function FGERIV. 


Function value: 

Not applicable. 


Labeled COMMON parameters: 
None. 


Blank COMMON parameters: 
None. 
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Inputs 
None . 

Outputs 

None. 

Storage Requirement 
Not applicable. 

Description 

FGERIV calculates the value of DERIV from the values of alpha and beta and 
the current radiance value. If DERIV is negative, it is reset to zero. 

If DERIV is zero, it is rese* to 1:E-10. 

Flowchart 

Not applicable. . . 

Listing 


See Appendix D for program. 


3.5.18 FUNCTION F6NCTN 
Purpose 

FGNCTN calculates the Y fit for the curve. 
Linkages 

FGNCTN Is called by CHNFIT. 

Interface 

Calling sequence: 

YFIT(I) = FGNCTN (X, I, A). 


Calling sequence parameters: 

Input/ 

Argument Output Descrl ption 


X 

I 

A 


I Array of time points. 

I Index to current point. 

I Array of coefficients A, Alpha, Beta and t^. 


Function value: 


Name Description 

YFIT(I) Evaluated function. 


Labeled COMMON parameters: 

Full descriptions of labeled COMMON blocks are contained In Appendix A. 

Element Input/ 

Label Variable Position Output Description 

/CHAN/ ACH 1 I A value returned from 

TRJFIT. 


/CHAN/ ALP 2 


Alpha value returned 
from TRJFIT. 


/CHAN/ ABET 3 


Beta value returned from 
TRJFIT. 
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Blank COMMON parameters: 

None. 

Inputs 

None. 

Outputs 

None. 

» 

Storage Requirement 
Not applicable. 

Description 

FGNCTN sets the alpha value positive and calculates F, current value of X/ 
alpha. If F Is less than or equa"' to 1., FGNCTN = ACH; otherwise FGNCTN 
Is calculated from the parameters of the training field. 

Flowchart 
Not applicable. 

Listing 

See Appendix D for program. 
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4.0 OPERATION 


4.1 OPERATING DESCRIPTION 

CLASFYT Is operational on the IBM 3031 computer at LARS, West Lafayette, 
Indiana. 

The CLASFYT program is one of the programs of the BADHWAR SYSTEM. 

CLASFYT requires the use of a D disk which i5 assigned as a temporary disk 
and an E disk which is used to hold LARS routines. The user therefore, 
must not assign a disk to his machine using either MODE E or MODE D. These 
disks will be assigned as needed. 

Prior to executing the CLASFYT program the user must establish on his A disk 
a USER INFORMATION file as described in Section 3.4.1. 
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4.2 COMMANDS DESCRIPTION 


To execute CLASFYT the user wfll enter the following series of commands 
which Invoke the JOB CONTROL SOFTWARE. These commands are divided Into two 
classes namely (1) FUNCTION coimiands and (2) PROGRAM commands . The FUNCTION 
commands, which perform all the functions except executing the program, are 
reusable, I.e., once they are Invoked they remain In effect until reentered. 
The PROGRAM commands, which simply execute the program, must be reentered 
each time the program Is to be executed. 

The following list gives the commands required to execute the CLASFYT 
program. They are all FUNCTION commands except the PROGRAM command CLASFYT. 
These commands are to be given In the listed order.. 

START 

CLASFYT 

END 

The following sections describe each of the commands In detail. Input 
fields are separated by blanks. If more than one word Is required to 
describe an Input field, the description Is enclosed In pointed brackets <>. 
If an Input Is optional the field Is enclosed In square brackets []. Do 
not Include these explanatory characters <> [] when actually submitting Input 
to the computer. To enter a command the user types one Input per defined 
Input field and separates each fie. with a blank. 
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4.2.1 START 


The START command spools the user's console file. The use of this comtand 
along with the END command will provide a listing of all Information 
appearing on the user console file. (If running an Interactive job this 
Is the terminal.) If running a batch job this Is a system defln-id device. 
The START command Is Invoked by the user typing the following: 

START 


4-3 


4.2.2 CLASFYT 


The CLASFYT conwand Is a PROGRAM command and is used to Invoke the execution 
of the CLASFYT program. All input Landsat data is from the LARS RT&E Data 
BASE and a series of programs is required to do the necessary interface. 
The following diagram illustrates this software flow. 

..RTEER <LARS R0UTINE> 

CLASFYT LCTINF ‘.SEGALL <LARS ROUTINE> 

:.TSRTHS....IVALUE 

For a detailed description of ttie JOB CONTROL SOFTWARE See Appendix C. 

This CLASFYT conmand must not be given unless the user has established a 
user's file on his A disc under the file description <SEGMENT NUMBER> USERT A 
as described in Section 3.4.1. Furthermore the SEGMENT NUMBER must contain 
4 digits of the segment. For example: segment 0882 dictates that the 

USER INFORMATION file be named: 

0882 USERT A 

The CLASFYT coimand is invoked by the user typing the following: 

CLASFYT <SEGMENT#> <ACQ.#1> <ACQ.#2> <ACQ.#3> <ACQ.#4> <ACQ.#5> 

As described in Section 3.4.2 the CLASFYT program requires the user to 
execute using 5 segment acquisitions. 

The classification file output from an execution of CLASFYT is written to 
a file named <SE6MENT NUMBER> CLASS A. Program inputs acquired from the 
LARS Data Base and the USER INFORMATION file is spooled to the HOUSTON 
printer. 


4.2.3 END 


This command closes the user's console file and causes a spooled copy to be 
sent to the HOUSTON printer. This command has no effect If the START 
command was not previously Issued. The END command Is Invoked by the user 
typing the following; 

END 


4.3 OPERATING EXAMPLE 


For our example we assume that the user has established on his A disk the 
required USER INFORMATION file. In addition since our example deals with an 
execution using segment 882 the user's Information file must be established 
under the filename 0882 and flletype USERT. 

COMMAND EXPLANATION OR ACTION TAKEN 

START Spools the console file. 

CLASFYT 0882 78159 78186 78222 78231 78267 

Executes the CLASFYT using 5 
acquisition of segment 882. 


END 


Closes the user's console 
file and prints the file. 










APPENDIX A 


COMMON BLOCK /LINES/ 


NXTLNE 

LINENO 

IDATA(856,8) 

NOACQS 


Not currently used. 
Current line number. 
Acquisition data. 

Number of acquisitions. 


A-1 


COMMON BLOCKS /MISC/ 


NOACQ 

ACQDT(5) 

KSEGM 

LANDST(5) 


Number of acquisitions. 

Acquisition dates stored left jus-^^-^ied in double precisi 
word. 

Segment number stored as 4 character word. 

Landsat numbers. 


COMMON BLOCK /USER/ 


ICROP(IO) 

ZQ(4,10) 

ZX(4,10) 

ZY(4,10) 

QZ(4,10) 

ISYMBL(IO) 

I0RDRC(4) 

NOCROP 


Names of crops of interest. 

Initial guesses for A, Alpha, Beta and t^^ for channel 

Initial guesses for A, Alpha, Beta and t^ for channel 

Initial guesses for A, Alpha, Beta and t^ for channel 

Initial guesses for A, Alpha, Beta and t^ for channel 

Symbols for crops of interest. 

Order of channel processing. 

Number of crops to process . 
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CLASFYT Variable Definitions 


ORIGINAL PAGE IS 
OF POOR QUALITY 

VARIABLE DEFINITION 

RCQDT BCOUISITIOM DRIES 
RLOG PERL LOG FUMCTTON 

RLPHl END CURVE FIT VRLUE CH 1» SET AFTER TRJFIT» USED BY BOTH CHNFIT 

RLPHS 2ND CURVE FIT VALUE CH S» SET AFTER TRJFIT» USED BY BOTH CHNFIT 

ALPH3 3RD CURVE FIT VALUE CH 3» SET AFTR TRJFITj USED BY BOTH CHNFIT 

ALPH4 2ND CURVE FIT VALUE CH 4f SET AFTER TRJFITj USED BY BOTH CHNFIT 

A1 INIT GUESS 1 CH 1j SET AFTER TRJFITj USED ALL CALLS TO CHNFIT 

R2 INIT GUESS 1 CH 2j SET AFTER TRJFITj USED ALL CALLS TO CHNFIT 

A3 INIT GUESS 1 CH 3j SET AFTER TRJFITj USED ALL CALLS TO CHNFIT 

A4 INIT GUESS 1 CH 4j SET AFTER TRJFITj USED ALL CALLS TO CHNFIT 

BETAi 3RD CURVE FIT VALUE CH Ij SET AFTER TRJFITj USED BY BOTH CHNFIT 

BETAS 3RD CURVE FIT VALUE CH £j SET AFTER TRJFITj USED BY BOTH CHNFIT 

BETAS 3RD CURVE FIT VALUE CH 3» SET AFTER TRJFITj USED BY BOTH CHNFIT 

BETA4 3RD CURVE FIT VALUE CH 4j SET AFTER TRJFITj USED BY BOTH CHNFIT 

ESEG SEGMENT NUMBER STORED LDGICAL^l 
CHNFIT CHANNEL FIT SUBROUTINE 
CHV HOLD AREA FDR CHANNEL VALUES 

CLNDST CONSTANTS FDR LANDSAT 

CPTIME SYSTEM TIME ROUTINE 

CRDPTM SUBROUTINE TO READ USER. INFORMATION FILE 
HBUF BUFFER FDR CLASSIFIED PIXELS 
HDR HEADER 
I TEMP. INDEX 

I COP RETURN PARAMETER FDR PIXEL U I THIN DEFINED VERTICES 
ICPTIM START TIME IN HUNDREDTH OF SECONDS 
ICPTIM INITIAL TIME 
I CROP NAME OF CROP 

ICUT MATRIX OF PIXELS CUTj INDEXED BY CHANNEL AND CROP 
ICI.ITT SUMM OF PIXELS CUT FDR EACH CROP 
I DATA BUFFER FDR RADIANCE VALUE TAPE 
IDfiY £ CHftRfiCTER DRY 

TDUPDT INDEX TO DUPLICATE DATE DR ZERO IF NO DUPLICATE DATE 
I HOLD 1 HOLD AREA FDR INTEGER 
IHDLD4 HOLD AREA FDR LOGICAL-^! CHARACTERS 
I MON 2 CHARACTER MONTH 

IDRDRC ORDER TO PROCESS CHANNELS FDR CHNFIT 

IPP PIXEL POINTER 

ISCREN NUMBER OF PIXELS WITH ZERO RADIANCE VALUES 
ISYM SYMBOLS FDR CROPS 

ITYPE SWITCH USED FDR CALL TO RDFLES FDR HERDER DR RAD. VALUES 
I YEAR £ CHAP CURRENT YEAR 
J TEMP. INDEX 

•JSEGM SEGMENT NUMBER AS DISPLAY "CODE 
K TEMP. INDEX 
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ORIGINAL bS 

VARIABLE DEFINITION QUALITY 

Kfl - KZ SYMBOL FDR P COMTIMUIHG MITH KB? KC» ... » KZ 
KBLPMK PLPHPMUMEPIC BLANK 
KCHRM CURRENT CHANNEL? USED IN DISPLAY 
KEOUAL ALPHANUMERIC EiIiUAL 

KSEGM SEGMENT NUMBER AS COMPUTATION NUMBER 
KUNCLS SYMBOL FDR UNCLASSIFIED PIXELS 
KURACQ INDEX TO CURRENT ACQUISITION 

KURCHN INDEX TO CURRENT CHANNEL 

KURCRP INDEX TO CURRENT CROP 

KIJRPIX INDEX TO CURRENT PIXEL IN TRAINING FLD 

KURPXL INDEX TO CURRENT PIXEL IN TEST FLD 

KZ THROUGH KZ AS THE SYMBOL FDR Z 

K1 - K4 SYMBOL FDR 1 CONTINUING WITH K£? K3? K4 
L TEMP. INDEX 

LACQDT ACQUISTIDN DATE AS LDGICAL^l 

LANDST ARRAY OF LANDSAT VALUES? DNE/ACQUISITIDN 

LIMIT MAXIMUM LIMIT FDR CALLS TO CHNFIT 

LINEND CURRENT LINE NUMBER FDR RADIANCE VALUES 

LLINE INDEX FDR LINE DO LOOP 

LDC TEMP LOCATION FDR UNPACKING INPUT DATA 

M TEMP. INDEX 

MAXLNE LARGEST LINE NUMBER IN CURRENT FLD 
MAXPIX LARGEST PIXEL NUMBER IN TRAINING FLD 
MAXPXL LARGEST PIXEL NUMBER IN TEST FLD 
MINLNE SMALLEST LINE NUMBER IN CURRENT FLD 
MINPIX SMALLEST PIXEL NUMBER IN TRAINING FLD 
MINPXL SMALLEST PIXEL NUMBER IN TEST FLD 
MODE WEIGHTING INDICATOR? ALWAY 1? NOT USED 
NFILES NUMBER OF ACQUISITION FILES 
NMAX MAXIMUM NUMBER OF CHOSEN PIXELS 
NDACQ NUMBER OF ACQUISITIONS 
NOACQS NUMBER OF ACQUISITIONS 
NDCRDP NUMBER OF CROPS TO CHECK 

NDKRDP NUMBER OF PIXELS THAT BELONG TO -NONE OF THE Cf'DPS 
NOP IX INDEX TO CURRENT PIXEL IN ORIGINAL PIXELS IN TRAIN FLD 
NPTSTF CURRENT NO OF POINTS IN TRAINING FLD 
HTERMS NUMBER OF DEGREES OF FREEDOM <4> 

NVCPTM ENDING TINE IN HUNDREDTHS OF SECONDS 
NVCPTM END INI? VIRTUAL TIME 
NXTLNE NOT USED 

PELTA STH ARGUMENT IN CALLING SEQUENCE TO TPJFIT? NOT USED OTHERWISE 

QCHI • CHI SQ VALUE RETURNED FROM TRAINING CALL TO TRJFIT FDR CH 1 
OCDEFS CH 1 1ST VALUE CURVE FIT APPROX? SET IN TRJFIT? USED BOTH CHNFIT 
QDAT CH 1 HOLD AREA FDR CHOSEN RADIANCE VALUES 

QFTX IITH PARAMETER IN CALL TO TRJPlT? NOT USED. 

QFTY IITH PARAMETER IN CALL TO TRJFIT? NOT USED. 

QISQR CH 1 CHI SQ VALUE FROM TRJFIT? USED IN BOTH CALLS TO CHNFIT 
QMAX SET EQUAL TO TCUT 

QMEAN CH 1 MEAN FDR DRIG. PIXELS? THEN SELECTED PIXELS 
QMEANN MEAN OF CHI SQ. VALUES AFTER TRJFIT AND CHNFIT FDR TRAINING FLD 

QPT SINGLE DIMENSION ARRAY FDR. CHAN 1 VALUES? USED TO CALL SDEV 

QR TEMP LOCATION 
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VARIABLE DEFINITION 


ORIGINAL PAC.:: 

OF POOR QUALITY 


?TFI DELTfi VftLUE IN CALL TO CHNFIT 

?TB DELTR VPLUE IN CFILL TO CHNFIT 

?DEXCT SUBROUTINE TO REBD' DflTFI FROM RTe<E DRTFI BFiSE 
5DFLES SUBROUTINE TO RERD RRDIRNCE VRLUES 

RECT SUBROUTINE TO DETERMINE IF PIXEL WITHIN DEFINED VERTICES 

SDEV SUBROUTINE TO CflLCULRTE MEFiN RND STRNDRRD DEVIRTIDN 

SGP SU3MP VALUE IN CALL TO CHNFIT 

SGB SIGMA VALUE IN CALL TO CHNFIT 

FOR ALL TEST PIXELS AND THEN CHOSEN PIXELS 

FOR ALL TEST PIXELS AND THEN CHOSEN PIXELS 

FOR ALL TEST PIXELS AND THEN CHOSEN PIXELS 

FOR ALL TEST PIXELS AND THEN CHOSEN PIXELS 

jMQj USED CALL TO TRJFIT» RECALC'D FDR TEST CHNFIT 

5MX» USED CALL TO TRUFITi. RECALC ^D FDR TEST CHNFIT 

;MYj 'USED CALL TO TRJFITs RECALC ^D FDR TEST CHNFIT 

;M2> used call to TRJFITi- PECALC^D FDR TEST CHNFIT 

SIGMA 9TH ARGUMENT IN CALLING SEQUENCE TO TRJFIT> NOT USED OTHERWISE 
SNO TEMP. VALUE IN CALCULATING CH 1 CONSTANT FDR CHNFIT 

SNX TEMP. VPLUE IN CALCULATING CH £ CONSTANT FDR CHNFIT 

SNY TEMP. VALUE IN CALCULATING CH 3 CONSTANT FDR CHNFIT 

SN2 TEMP. 

SO. ST. DEV 

SORT SQUARE ROOT FUNCTION 

ST TEMP LDC ST. DEV. FDR TRAIN FLD APPROX AFTER CHNFIT CALL CHAN> 

SX ST. DEV. FDR CH £ FDR CHOSEN PIXELS <SGMX) j THEN CHI SQ<NDT USED) 

SY ST. DEV. FDR CH 3 FDR CHOSEN PIXELS fSGMY) » THEN CHI SQ<NDT USED) 

S2 ST. DEV. FDR CH 4 FDR CHOSEN PIXELS <SGMZ) j THEN CHI SQ^NDT USED) 


SGMQ 

CH 

1 

ST. 

DEV 

SGMX 

CH 

o 

L. 

ST. 

DEV 

SGMY 

CH 

3 

ST. 

DEV 

SGM2 

CH 

4 

ST. 

DEV 

SGMl 

CH 

1 

FTN 

OF 

SGM2 

CH 

£ 

FTN 

OF 

SGM3 

CH 

J 

FTN 

OF 

SGM4 

CH 

4 

FTN 

OF 


VALUE IN CALCULATING CH 4 CONSTANT FDR CHNFIT 

FDR CH 1 FDR CHOSEN PIXELS <SGMO) » THEN CHI SQ<NDT USED) 


TOUT CONSTANT 6. £3? RESET TO 7.036 IF DUPLICATE DATE 

TFTOIM CHI MEAN TRAIN FLD PIXEL TO'S <TRPT01) AFTER CHNFIT ON TRAIN FLD 
TFT0£M CH£ MEAN TRAIN FLD PIXEL T0"S <TRPT0£) AFTER CHNFIT ON TRAIN FLD 
TFT03M CH3 MEAN TRAIN FLD PIXEL TO^S <TPPT03) AFTER CHNFIT ON TRAIN FLD 
TFT04M CH4 MEAN TRAIN FLD PIXEL TO'S <TRPT04) AFTER CHNFIT ON TRAIN FLD 
TODAY TODAY ^S DATE 


TRAJECTORY FIT SUBROUTINE 

TO VALUE FDR CH 1 FDR EACH PIXEL AFTER CHNFIT ON TRAIN FLD 

TO VALUE FDR CH £ FDR EACH PIXEL AFTER CHNFIT ON TRAIN FLD 

TO VALUE FDR CH 3 FDR EACH PIXEL AFTER CHNFIT ON TRAIN FLD 

TO VALUE FDR CH 4 FDR EACH PIXEL AFTER' CHNFIT ON TRAIN FLD 

LOWEST VALUE FDR TEST PIXEL TO <TPAIN FLD TO MEAN +-.£5 DR .£ 
HIGHEST VALUE FDR TEST PIXEL TO (TRAIN FLD TO MEAN +-.£5 DR .£ 


TRJFIT 
TRPTOl 
TRPT0£ 

TRPT03 
TRPT04 
TRTOLW 
TRTOUP 

TRTOIM TRAINING FIELD TO MEAN FDR CH 
TRT0£M TRAINING FIELD TO MEAN FDR CH 
TRT03M TRAINING FIELD TO MEAN FDR CH 
TRT04M TRAINING FIELD TO MEAN FDR CH 
TOl TO VALUE RETURNED FROM TRJFIT 
T0£ TO VALUE RETURNED FROM TRJFIT 
T03 TO VALUE RETURNED FROM TRJFIT 
T04 TO VALUE RETURNED FROM TRJFIT 
T1 FTN ST. DEV CHOSEN PXLS<SGMQ) 
TlX TEMP USED IN CALCULATING SGMl 
T2 FTN ST. DEV CHOSEN PXLS <SGMX) 
T£X TEMP USED Iii CALCULATING SGM£ 


1 

£ 

O 

4 

FDR CH 1 
FDR CH £ 

FDR CH 3 
FDR CH 4 

DR CURRENT PXLj USED TO CALC SGMl 
FDR BOTH CHNFIT CALLS 
DR CURRENT PXL? USED TO CALC SGM£ 
FDR BOTH CHNFIT CALLS 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


VARIABLE DEFINITION 

T3 FTN ST. DEV CHOSEN PXL3<SGMY> OR CURRENT PXL» 

T3X TENP USED IN CPLCULBTITlG SENS FOR BOTH CHNFIT 

T4 FTN ST. DEV CHOSEN PKLS <Si3NZ> OR CURRENT PXL> 

T4X TEMP USED IN CBLCULBTING SGM4 FOR BOTH CHNFIT 

HEIGHT I.IEIGHTING VALUES l./SIGMfl^^a> 

XCH HOLD BREA USED IN UNPACKING RADIANCE VALUES 
XCHI CHI SI? VALUE RETURNED FROM TRAINING CALL TO TRJFIT FOR CH 2 
XCHIN CONTROL VALUE FOR CALLS TO TRJFIT AND CHNFIT 

XCOEFS CH £ 1 ST VALUE CURVE FIT APPROX y SET IN TRJFIT, USED BOTH CHNFr 

XDAT CH £ HOLD AREA FOR CHOSEN PIXEL 

XDAV ARRAY OF ACOUISITION DAYS OF YEAR / 100. 

XDIF TEMP 

XDIFTM DIFFERENCE IN STARTING AND ENDING TIME 
XDIFVT DIFFERENCE IN STARTIN'o AND ENDING VIRTUAL TIMES 
XFD USER SPECIFIED X COORD I ANTES OF FIELDS 

XFIT IITH PARAMETER IN CHNFIT, NOT USED OTHERWISE 

XHl CHI HIGH LIMIT FOR TRAIN PIXELS, MEAN OF DRIG PIXELS + 3 ST DEV 

XH£ CHS HIGH LIMIT FOR TRAIN PIXELS, MEAN OF ORIG PIXELS + 3 ST DEV 

XH3 CHS HIGH LIMIT FOR TRAIN PIXELS, MEAN OF ORIG PIXELS + 3 ST DEV 

XH4 CH4 HIGH LIMIT FOR TRAIN PIXELS, MEAN OF ORIG PIXELS + 3 ST DEV 

XISiT'R CH £ CHI SQ. VALUE FROM TRJFIT, USED IN BOTH CALLS TO CHNFIT 

XLPMDA TOLERANCE VALUE USED IN CALLS TO TRJFIT AND CHNFIT 

XLl CHI LOW LIMIT FOP TRAIN PIXELS, MEAN OF ORIG PIXELS - 3 ST DEV 

XL£ CHS LOW LIMIT FDR TRAIN PIXELS, MEAN OF DRIG PIXELS - 3 ST DEV 

XL3 CHS LOW LIMIT FDR TRAIN PIXELS, MEAN OF DRIG PIXELS - 3 ST DEV 

XL4 CH4 LOW LIMIT FDR TRAIN PIXELS, MEAN OF DRIG PIXELS - 3 ST DEV 

XMAX SET EQUAL TO TCUT 

XMEAN CH £ MEAN FDR DRIG. PIXELS, THEN SELECTED PIXELS 
XMEANN MEAN OF CHI SQ. VALUES AFTER TRJFIT AND CHNFIT FDR TRAINING FLD 
XPT SINGLE* DIM ARRAY FDR CH £ RADIANCE VALUES, USED TO CALL SDEV 
XR TEMP 

XTJ TEMP HOLD AREA FDR FLOAT PT DATE VALUE 
XTDL MINIMUM TOLERANCE AFTER CALL TO CHNFIT 

VCHI CHI SQ VALUE RETURNED FROM TRAINING CALL TO TRJFIT FDR CH 3 

YCDEFS • CH 3 1ST. VALUE CURVE FIT APPROX, SET IN TRJFIT, USED BOTH CHNFIT 

YD AT CH 3 HOLD AREA FDR CHOSEN PIXEL 

YFD USER SPECIFIED Y COORDINATES OF FIELDS 

YFT IITH PARAMETER IN CALLS TO CHNFIT, NOT USED OTHERWISE 

YISQR CH 3 CHI SQ. VALUE FROM TRJFIT, USED IN BOTH CALLS TO CHNFIT 

YMAX SET EQUAL TO TCUT 

YMEAN CH 3 MEAN FDR DRIG. PIXELS, THEN SELECTED PIXELS 
YMEANN MEAN OF CHI SQ. VALUES AFTER TR-JFIT AND CHNFIT FDR TRAINING FLD 

YPT SINGLE DIM ARRAY FDR CH 3 RADIANCE VALUES, USED TO CALL SDEV 

YR TEMP 

2A TEMP APPROX FIT VALUE CHNFIT IN TRAIN FLD, RESET IN CHNFIT TEST 

ZB TEMP APPROX FIT VALUE CHNFIT IN TRAIN FLD, RESET IN CHNFIT TEST 

2CHI CHI SQ VALUE RETURNED FROM TRAINING CALL TO TRJFIT FDR CH 4 
2CDEFS CH 4 1ST VALUE CURVE FIT APPROX, SET IN TRJFIT, USED BOTH CHNFIT 

2D AT CH 4 HOLD AREA FDR CHOSEN PIXEL 

ZFL LINE NUMBER USED IN CALL TO RECT 

2FP PIXEL NUMBER USED IN CALL TO RECT 


USED TO CALC SGM3 
CBLLS 

USEd'^TD calc SGM4 
CALLS 
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VARIABLE DEFINITION 

I COP CH 4 CHI SO. VHLUE FROM TPJFITi USED IM BOTH CHLLS TO CHNFIT 
MFt^< SET EQUAL TO TCUT 

MEAN CH 4 MEAN FDR OPIG. PIXELS^ THEM SELECTED PIXELS 
MEAMM ME AM OF CHI SO. VALUES AFTER TRJFIT AMD CHMFIT FOR TRAIMING FLD 
PT SIMGLE DIM ARRAY FOP CH 4 RADIAMCE VALUES? USED TO CALL SDEV 
0 IMITIAL USER GUESS FDR CH I CURVE 

P TEMP LDCATIDM 

X IMITIAL USER GUESS FOP CH 2 CURVE 

V IMITIAL USER GUESS FDR CH 3 CURVE 

Z IMITAIL USER GUESS FDR CH 4 CURVE 
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